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Abstract 

We propose an extension of Dubrovin's perturbative approach to the study 
of normal forms for non-Hamiltonian integrable scalar conservation laws. The 
explicit computation of the first few corrections leads to the conjecture that such 
normal forms are parametrized by one single functional parameter, named vis- 
cous central invariant. A constant valued viscous central invariant corresponds 
to the well-known Burgers hierarchy. The case of a linear viscous central in- 
variant provides a viscous analog of the Camassa-Holm equation, that formerly 
appeared as a reduction of a two-component generalization of the Camassa- 
Holm equation. We write explicitly the negative and positive hierarchy asso- 
ciated with this equation and prove the integrability showing that they can 
be mapped respectively into the heat hierarchy and its negative counterpart, 
named Klein-Gordon hierarchy. A local well-posedness theorem for periodic 
initial data is also proven. Finally, we explicitly construct asymptotic solutions 
obtained by quasi-Miura transformations and discuss Dubrovin's Universality 
conjecture for the local behaviour of the solutions near the point of gradient 
catastrophe. 

1 Introduction 

In [12] B. Dubrovin proposed a perturbative approach to the study of integrable 
Hamiltonian evolutionary equations (or systems of equations) of the form 

ut = A{u)ur, + e[Bi{u)uxx + B2{u)ul] + e^[Ci{u)uxxx + C2{u)uxUxx + C:i{u)ul] + C(e^). 

(1.1) 
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where e is a small formal expansion parameter. If the unperturbed equation obtained 
for e = is integrable and Hamiltonian (this is always true in the scalar case), 
Dubrovin's method provides an effective tool to construct perturbative dispersive 



equations of the form (1.1) preserving integrability and Hamiltonian properties. The 



r.h.s. of (1.1 ) is a formal expansion, it contains in principle infinitely many terms and 



no additional assumptions concerning convergence are usually enforced. It is natural 



to classify equations of the form (1.1) up to the group of Miura transformations 
defined via formal series of the form 



V = f{u) + t[g{u)uj^ + t^[hi{u)ux^ + h2{u)ul] + ... 



:i.2i 



Since most of physical applications involving equations of the form (1.1) contain a 



finite number of terms, one would be tempted to work with truncated Miura trans- 
formations only. This choice turns out to be ineffective and the reason is two-fold. 
First of all, an infinite formal expansion is needed, in general, to define the inverse 
of a Miura transformation, truncated or not: if one does not allow infinite formal 
expansion, then the set of truncated Miura transformations do not form a group. 
Secondly, the restriction to finite expansions will exclude from the analysis all those 



non-evolutionary equations that can be put in the evolutionary form (1.1). Let us 
consider for example the celebrated Camassa-Holm equation 



ut- e Uxxt 



-?)UUx + e [uu^ 



;i.3) 



Introducing the transformation 



u 



its formal inversion is given by an infinite formal series that allows to re-write the 



equation (1.3) in the evolutionary form (1.1) with respect to the variable v 



However in a number of cases truncated equations might have their own interest 
and in particular those equations possessing infinitely many approximate symmetries 



In this paper we will extend Dubrovin's approach to the case of integrable conserva- 
tion laws of the form 

Ut = dx {u^ + e[a{u)u4 + e^[bi{u)uxx + b2{u)ux2] + C(e^)} . (1.4) 

Equations of this form are not necessarily Hamiltonian. For instance, in the scalar 
case, the operator dx defines a Poisson structure on the space of functionals F[u] = 
J f{u,Ux, Uxx, ■ . .) dx, but in general the current 

+ e[a{u)ux] + e^[bi{u)uxx + &2(m)mx2] + ^(e^) 
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is not the variational derivative of a functional F[u\. To overcome this difficulty we 
consider the extension of the Hamiltonian formalism to 1-forms introduced in [1]. 
Within this framework, currents are viewed as 1-forms in a suitable jet space on 
which the operator dx defines a Poisson structure. Hence, the extension of Dubrovin 
perturbative scheme is immediate apart from the choice of the class of admissible 
formal Miura transformations. There are two natural options: 

1. To consider canonical transformations as in jl2j only . This choice is motivated 
by results in [15] that guarantees any deformation of the operator to be 
eliminated by a Miura transformation. 



2. To consider Miura transformations that preserve the form of the equation ( 1.4). 
This class is clearly wider and plays a crucial role in the elimination of few 
inessential functional parameters. This is the reason why we have chosen to 
work with this second option. 

An alternative approach to the classification of integrable hierarchies, based on the 
nonlinear perturbations to linear equations, has been developed for instance in [27], 
[2H], [23, [2]- This is the so called symmetry approach. 

The main results of this paper can be summarized as follows: 

• We classify up to the fifth order in the deformation parameter scalar integrable 



conservation laws of the form (1.4). We will focus on the viscous case corre- 
sponding to a{u) 7^ and for which the corresponding currents are not exact 
1-forms. It turns out that, up to the class Miura transformations preserving 



the form of the equation (1.4), all deformations are uniquely determined by 
the function a{u) named viscous central invariant. For this reason, we refer to 
the corresponding equations as integrable viscous conservation laws. Dispersive 
deformations obtained for a{u) = and bi{u) ^ have been already studied in 



in the case of Hamiltonian conservation laws . 
We observe that except for the case of constant viscous central invariant where 



the right hand side of (1.4) becomes finite and gives the well known Burgers 
equation, the generic integrable conservation law contains infinitely many terms. 
We show that the case of a linear viscous central invariant a[u) = u corresponds 
to the equation 

Ut - euxt = dx {u^ - euux) , (1.5) 

formerly derived as a non-Hamiltonian reduction of two-component generaliza- 
tions of the Camassa-Holm equation [3T1 [16] . A few interesting properties of 



equation (1.5) have been already investigated, such as the existence of a typ- 



ical viscous shock dynamics associated with cliffons, i.e. weak traveling wave. 
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solutions [16]. Nevertheless, its integrable structure has not been yet studied 



in detail. We observe that, remarkably, equation (1.5) is related by a Miura 
transformation to the equation 



Ut 

which is a symmetry of the equation 



a. ^X^e^MM(fc)j (1.6) 



Ur^,=dA- + ed^^] (1.7) 



u u 



appearing in Calogero's list [4J. We show that the equation (1.7) is the first 
member of the full integrable "negative" hiearchy 

Ur^„ = dJ- + ed,-\ (1) n = 1,2,3,... (1.8) 



u u ^ 

We also construct the positive counterpart which turns out to be non-local, that 



is its evolutionary form is not truncated, as it happens for the equation (1.6) 



It should be noted that the negative hierarchy (1.8) and its positive counter- 
parts can be obtained as a one component reduction on the negative hierarchy 
recursively constructed in |16]. We provide an alternative proof of the integra- 
bility by showing that the negative and positive hierarchies can be transformed 
respectively into the Burgers and then the heat hierarchy and its negative coun- 
terpart. Interestingly the first member of the negative counterpart is given by 
the Klein-Gordon equation and for this reason we refer to it as Klein-Gordon 
hierarchy. 

We study the well-posedness (existence, uniqueness and continuous dependence 



on initial data) for the periodic Cauchy problem for the equation (1.5) (with 
e = 1) locally in time, using the general techniques developed by Kato in [18], 
similarly to what has been done for the Camassa-Holm equation in [7]. More 



specifically, we prove that the periodic Cauchy problem for the equation (1.5) 
is locally well-posed for m G , where m = u — u^- It is a work in progress 
the determination of the conditions that guarantee the global well-posedness of 
the periodic Cauchy problem and the blow-up conditions in finite time, in full 
analogy with the results obtained for the Camassa-Holm equation [71 [6] . We 
note that in [12] well-posedness and blow-up properties have been studied for 



the two component generalization of the equation (1.5) introduced in [5] and 

HSl. 



Using Mathematica [3l] we investigate numerically the Cauchy problem for the 



equation(1.5) assuming periodic boundary conditions. This has been done by 



4 



converting the equation (1.5) into a system consisting of an evolutionary PDE 
together with an ODE. The numerical integrator, based on the method of lines 
for solving the PDEs, applies a pseudo-spectral method or a spatial discretiza- 
tion with finite differences of the 4th order. We perform numerical simulations 
with very simple initial data such as translated sinusoidal waves with different 
wave lengths and observe that solutions behaviour is as expected: the initial 
datum evolves in low viscosity regime approaching the gradient catastrophe 
and the steep profile so created is asymptotically dissipated. We also provide 
numerical evidence of the existence of initial conditions for which the solution 
breaks in finite time similarly to what happens for the Camassa-Holm equation 
in [7]. 

• The class of perturbed equations considered in this paper coincides up to O(e^) 
with the generalized Burgers equations studied in [9]. We show that near the 
critical point the solutions of scalar Hamiltonian conservation laws should sat- 
isfy a second order ODE which is the non-Hamiltonian analogue of the Painleve 
equation arising in the description of the critical behaviour of Hamiltonian con- 
servation laws [12]. The mentioned ODE admits the Pearcey function as a 
particular solution. This is consistent with the results in [21] and [9] concerning 
the universality of the critical behaviour of Burgers' equation and its general- 
ization. 

The paper is organized as follows. In Section 2 we give the definition of integrable 
scalar conservation laws and we introduce their normal form in Section 3. In Sec- 
tion 4 we present the perturbative classification results concerning integrable viscous 
conservation laws and discuss in detail the case of a linear viscous central invariant 
in Section 5 . Section 6 is devoted to the proof of the local well-posedness for the 
Cauchy problem and numerical simulations are presented in Section 7. In Section 8 
we review some results about the quasitriviality of evolutionary PDEs due to Liu and 
Zhang [23]. However, instead of characterizing the quasitriviality transformations in 
terms of their infinitesimal generators as in [23], we provide a direct characterization 
of the coefficients of the transformation. In Section 9 these results are applied to the 
study of the critical behaviour of solutions near the critical point. 



2 Commuting flows 

In the present section we introduce the main notions concerning integrability and 
perturbation theory of conservation laws. For the sake of simplicity we specify all 
definitions and properties in the case of scalar conservation laws, that is the subject 
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matter of this paper. The interested reader will find a rigorous and more general 
treatment in the paper [1]. 
Let 

ut = d^a{u,u^,Uxcc, ■ ■ ■) (2.1) 

be a scalar conservation law. If we assume that the current a is a differential poly- 
nomial, then, rescaling the variables as a; — )■ ex,t — et, we can write the right hand 
side of any scalar conservation law as 

ut = ^o{u) + ^ e^ak{u, m^, u^^, ...)^ (2.2) 

where are differential polynomyals of degree k. Here the degree is assigned ac- 
cording to the following rule: deg/(M) = 0, degU(fc) = k. Clearly if we start from a 
differential polynomial, N is finite. However as mentioned above we also include the 



case = oo. In this case the equation (2.2) will be referred to as formal conservation 
law. 



Definition 2.1 A (formal) scalar conservation law (2.2) is said to he integrable if it 
admits infinitely many symmetries of the form 

Ur = d^(3{u,u^,u^^,...). (2.3) 

A very important subcase is given by the class Hamiltonian conservation laws. In 
this case the current a is the variational derivative of a suitable local functional (the 
Hamiltonian functional) 

H[u] = J h{u,U:i:, ...) dx 

where the function h is the Hamiltonian density. In this case we have 



^_5H _dh _^ / dh\ g2 ( dh 
5u du ^ \dux) ^ \duxx 

Two functionals F[u\ and G[u\ such that {F, G} = are said to be in involution. 

It is well-known (see e.g.[10J) that the operator dx defines a Poisson bracket in the 
space of local functionals. Given two local functionals F[u] and G[u], their Poisson 
bracket is the local functional defined by 

{FM.GH} = /|^a.(^).x. (2.4) 

It was proved in [191 E] that any local Hamiltonian operator can be reduced in the 
scalar case to dx by means of a Miura transformation. This means that any scalar 
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Hamiltonian equation (w.r.t a local Hamiltonian operator) can be written as a conser- 
vation law after a suitable Miura trasformation. Moreover, in the Hamiltonian case 
also the symmetries of the equation are Hamiltonian and the corresponding Hamil- 
tonian functionals are in involution with the Hamiltonian functional of the original 
equation. 

The previous analysis shows that the function a in (2.2) can be thought as a 
1-form. Such a 1-form is in general non-exact and the equation is not Hamiltonian 
with respect to the Hamiltonian operator d^- 

Nevertheless, extending a well-known construction in the finite dimensional set- 
ting, it is possible to define a Poisson bracket on the space of 1-forms Ai. In the 
scalar case, this is given by: 



da 



d(3 
dui 



0. 



(2.5) 



The above bracket satisfies the following properties (see [I] for more details): 

1. If a = 5F and /3 = 5G, then {a, /3} = G}. 

2. {•, ■} equips the space of 1-forms Ai with a Lie algebra structure; 

3. the Poisson structure induces an (anti)-homomorphism of Lie algebras between 
(Ai, {■,■}) and the space of evolutionary vector fields equipped with the Lie 
bracket given by the Lie commutator. 

From the last property it follows that two scalar conservation laws commute if and 
only if the associated currents are in involution. For convenience of the reader we 
give an alternative elementary proof of this fact. 



Proposition 2.2 Two flows of conservation laws of the form (2.2) and (2.3) com- 
mute iff a and /3 are in involution with respect to the bracket (2.5). 

Proof: Let us introduce the function if such that u = (fx- Then, integrating once 
w.r.t. the variable x, equations (2.2) and (2.3) give 

Lpt =a{u, Ux, Uxx, ■■■) + fit, t) 

ifr =f3{u, Ux, Uxx, • ■ • ) + git, ^) 

where / and g are two arbitrary functions. The request that the two equations above 
commute for any solution ip, i.e. 
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implies that the two compatibihty conditions 

ttr = Pt fr = 9t 

must hold separately. The second condition is identically satisfied if one choses f = ipt 



and g = ipr- Using the chain rule together with the equations (2.2) and (2.3) the first 
condition can be equivalently written as follows 



'-^ OU(i) out 



J 

The proposition is proved. ■ 
As a consequence of the above proposition, we have that any integrable hierarchy 

of conservation laws is always defined by an infinite family of 1-forms in involution. 

From this perspective Hamiltonian and non-Hamiltonian conservation laws can be 

understood within the same framework. 

Example Let us consider the Burgers hierarchy defined as 

ut„ = d^Un = [{u + d^Y u] , n = 0, 1, 2, ... (2.6) 



The fiows associated with ujn are in involution w.r.t. the Poisson bracket (2.5), i.e. 

{Wn, Wm} = 0. 

3 The normal form of conservation laws 

Let us consider a formal conservation law of the form 

ut = {u^ + ea{u)u^ + e^{bi{u)u^^ + b2{u)ul) + 0(e^)) . (3.1) 

It can be viewed as a higher order perturbation of the Hopf equation 

Ut = dxU^ = 2uUx- (3.2) 

Let us observe that there is no loss of generality by taking the Hopf equation as the 
leading order. Indeed, any equation of the form Ut = K{u)ux can be transformed into 



the equation (3.2) by a re-parametrization of the dependent variable u = u{v). The 



Hopf equation (3.2) is Hamiltonian in the standard sense and completely integrable. 
Hence, there exist an infinite set of symmetries parametrized by a function of one 
variable f{u) of the form 

Ut = d^f{u) = f'{u)u^. (3.3) 

It is straightforward to check that dfUr = drUt for any function f{u). 

The aim of this section is to illustrate a perturbative approach to the problem of 



classifying integrable deformations of the Hopf hierarchy of the form (3.1). Various 
approaches have been developed to construct such deformations, and they can be 
briefiy outlined as follows: 
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• Symmetries approach: It is based on the requirement that symmetry of the hy- 
drodynamic hmit can be extended to all perturbative orders |2l|33j. This is the 
most general approach because it does not require any additional assumption on 
the hierarchy. Nevertheless, it turns out to be computationally very expensive, 
in general. 

• Hamiltonian approach: proposed by Dubrovin in [12] it consists in deforming 
Hamiltonian densities of the Hopf hierarchy with the request that deformed 
Hamiltonian functionals are still in involution. 

• bi- Hamiltonian approach: proposed by Dubrovin and Zhang in [13] for the study 
of integrable PDEs appearing in the theory of Gromov- Witten invariants [2U [23] 
it is based on the extension of the bi-Hamiltonian structure to all perturbative 
orders. 

The approach we propose here is a generalization of the Hamiltonian approach. 
Instead of deforming the Hamiltonian functionals we deform their differentials. More- 
over, no assumption on the exactness of the deformed 1-forms is necessary. In this 
way, as mentioned in the previous section, we are able to treat simultaneously Hamil- 
tonian and not Hamiltonian integrable deformations. Since the former have been 
already extensively studied in [12], we will mainly focus on the latter. 

The classification procedure discussed in this paper is based on the following: 

Definition 3.1 ^ conservation law associated with the 1-form 

oo 

<^f^ = fiu) + ^ (^""gniu, U,^, ...) 
n=l 

where each gj is a homogenous differential polynomial of order j , is said to be inte- 
grable up to the order if there exists a 1-form 

oo 

i^f^ = fiu) + ^ e"^„(M, Mx, • • • ) 

n=l 

such that their Poisson bracket {u'l^^ .oj'f^} = vanishes identically modulo 0{e'^^^) 
for any function f{u). 

Let us observe that taking for instance f{u) = m^, the request that the two forms 
io'^t^ and w^^'^ commutes for any f{u) up to a certain order in e means that the 
corresponding deformed conservation law inherits all the symmetries of the Hopf 
equation up to the same order. 
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In the case of exact 1-forms the above commutativity condition reduces to the 
standard commutativity condition between Hamiltonian functionals with respect to 
a canonical Poisson bracket (see e.g. 



A direct apphcation of the definition (3.1 ) provides an effective tool to classify ap- 
proximate integrable conservation laws viewed as deformations to the Hopf equation. 
The explicit derivation of integrability conditions becomes computationally more and 
more expensive as the order of such deformations increases. Nevertheless, it turns out 
that the most general deformation of the Hopf equation contains a certain number of 
redundant functional parameters that can be eliminated using the invariance of the 
form of conservation laws with respect to a special class of Miura transformations. 

For instance, let us consider a general deformed conservation law of the form 



Ut 



(3.4) 



with 



i^'^f = u^ + ea{u)u^ + (6i(m)u^.^ + h2{u)ul) + {ci{u)u^^^ + C2{u)u^u^^ + C3(u)m^) 
+ ((ii(u)M4x + d2{u)u:r.Uxxx + dj,{u)ul^ + di{u)ulun,x + d5{u)ul) + (ei(u)u5a; 
+e2{u)uxUix + e^{u)uxxUxxx + eA{u)uluxxx + e^{u)uxul^ + ee{u)uluxx + e7{u)ul) 
+ 0{e% 



We observe that the form of the equation (3.4) is preserved under the Miura trans- 
formation 

u ^ V = u + e^dx(3{u^ Ux, ■ ■ ■) (3.5) 

where /3 is a homogeneous differential polynomial of degree k — 1. Hence, applying 
the Miura transformation (3.5) to the equation (3.4) we get 

Vt = dx {u^2{u{v,Vx, ...),Uxiv,Vx, ...), ...) + (3t{u{v,Vx, ...),Uxiv,Vx, ...), ...)) , 
where 



u{v,Vx, ...) 



e''dx(3iv,Vx,...) + 0ie''+'). 



Clearly the equation truncated at the order k — 1 coincides with the original equation 
truncated at the same order. At the order k we have 



fc-i 



ujk-^ujk = uJkiv,Vx, ...) - 2vdx(3{v,Vx, ...) + 



d/3{v,Vx,...)Q,+^^_2 



or equivalently 



k~l 

ujk{v,vx, ...) + 



dP{v,Vx, ...] 



E 

.1=1 



s=0 



s + 1 

I 



dv 



is) 



V[l)V(^s+l-l) 



(3.6) 
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For instance, for /c = 2 it is /3 = /3iVx and 



= hv^x + (62 + 2/3i)t;^ 

Hence the term v"^ can be eliminated by choosing /3i 

/3 = /32iVxx + and 



-62/2. For /c = 3 we have 



f^s = ci(t;)t;^^^ + (c2 + Q(^2i)vxVxx + (ca + 4/322)^^: 



Analogously, both terms VxVxx and can be eliminated. Similarly, it is straight- 
forward to check that for k = A one can retain just the terms V/^x and VxVxx) or 
alternatively f^^, while for = 5 only v^x and VxVix, or alternatively Vxx^xxx^ will 
survive after a suitable Miura transformation. 

The above examples suggest that it should be possible to choose [3 in such a way 
that only those terms that do not factor through Vx survive after the action of a 
suitable Miura transformation. This turns out to be true: 



Theorem 3.2 Given the conservation law 

Ut = dxUu^ 



(3.7) 



with = + ea{u)ux + 'Yl,k>2^^^k{u,Ux, ■ ■ ■) there exists a sequence of Miura 
transformations of the form 

u V = u + e'^dx(3{u, Ux, ■ ■ ■) 



that brings equation (3.7) into the normal form 

vt = dxUJ^2 

where 



{3.t 



UJ,,2 = V 



+ ea{v)vx + ^ e^ujk{v, Vx, 

k>l 



is such that 



dujk_ 
dvx 



0, VA; > 1. 



Proof: Let us first introduce some useful notations. A general 1-form /3 homoge- 
neous of degree k — 1 will be denoted as follows: 

P[il,i2,-M-i\''^{l) ■ ■ ■ "^(fc-l)- 

h,...,ik-i ii+2i2-\ h(A;-l)ifc_i=A;-l 

We will prove that the coefficients /3[ji,i2,...,ifc_i] can be are recursively determined 
through a lower triangular relation with respect to a suitable ordering of the mono- 
mials f . . . ti^^^i)- We introduce such an ordering, which is similar to a reverse 
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lexicografic ordering: the monomial f ^^-^ . . . ranks higher (or comes first or has 



(k-i) 



a higher rank) than v^|l^ . . . symbohcally 



(fe-i)' 



^(fc-1) ^ ^(1) • • • ^(fc-l)' 



if there exists m G 1} such that ii = ji for alH > m and im > jm- In other 

words, monomials are ranked according to the higher derivative of v. For instance, 
among the homogenous monomials of degree k — 1, the highest ranking monomial 
is f(/c-i) and the lowest ranking is v'^i^^- This gives clearly a total ordering of the 
monomials. In the following, we will write down a 1-form homogeneous of a certain 
degree from the highest to the lowest ranking. For instance, for a homogeneous 1-form 
of degree 3 

13 = /3[0,0,lF(3) + I3[1,1,0]V{1)V(2) + l3[3fl,0]vfi)- 

Let us now fix k and write Wfc as a + f (i)?7 where t] is ordered as above. In the 
rest of the proof a will be neglected as we will focus only on the terms of f (i)^] that 
factors through t>(i). 

Let us now expand the expression 



fc-i 

E 

s=l 



df3{v,v^,...) 



E 

.1=1 



I 



(3.9) 



entering the Miura transformation (3.6) as follows 

k-l 



E 



k-l 



dv 



s=3 



s-l 

E 

.1=2 



S+1 
I 



V(l)V(s+l-l) 



where 



fc-i 



^ = 2^ 2(s + l)v(s) 



s=i - '^.^) 

is meant to be ordered as above and the residual term is 



fc-i 



*=E 



s=3 



's-l 

E 

.1=2 



V(l)V(s+l-l) 



(3.10) 



Notice that the coefficients in /3 coincide, up to a positive factor, with the coefficients 
in (3, that is 



*2, ■■•,«*-!] 



where Cj^_..,,jj._j are positive constants. In particular, terms in /3 have the same ranking 
of the corresponding terms in /3. 
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We are now able to prove that it is always possible to eliminate all terms containing 
f (1) starting from the highest ranking coefficients. However, special attention should 
be paid as one needs to take into account contributions coming from both /3f (i) and 
rjvii) and 5. It will turn out that all such contributions appear in a lower triangular 
form and then we can recursively fix coefficients of /3 and consequently /3 in such a 
way to kill ?7f(i). 

Let us first note that is is always possible to fix the highest ranking coefficient 
in /3f (1) that annihilates the highest ranking coefficient in f (i)?]. Indeed, the highest 
ranking term in /3f (i) is /3[o,o,...,i]'^(i)'^(fc-i) while the highest ranking term in 77^(1) is 
'7[o,o,...,i]'y(i)^{fc-i)- It is immediate to check that no such a term can originate from 5 
and therefore /3[o,o,...,i] is uniquely determined by the condition /3[o,o,...,i] +^[o,o,...,i] = 0. 

Let us now proceed by induction following the ordering (from the highest to the 
lowest ranking) of the monomials. Assuming that the L highest ranking monomial in 
/3 is determined, we show that it is possible to find the L + 1 highest ranking monomial 
■ ■ ■ "^(fc-i) ^^^^ annihilates the corresponding L + 1 highest ranking 
monomial V[h,--.,ik-i\'^''(\)^ ■ ■ -"^(fc-i) by imposing the condition 

^[n,...,ifc_i]^^;}^^ • • • + hh,-,ik-i]'"'{iT ■ ■ ■ ^(fc-i) + residual terms = 0, (3.11) 

where residual terms come from 5. Indices ii, . . . , ik-i above are fixed. 

From its definition, 5 can be thought as an operator acting linearly on /3, that is 
5{(3i + [32) = S{f3i) + S{f32). Therefore, without loss of generality, we can focus our 
attention on the action of 6 on the generic term of the form 

Up to a positive factor, we can write 

5(7) = E E • • • • • • • • • -sr' • • • ^S-J)- (3.12) 

s=3 1=2 

Since I < s — 1 < s and analogously s + 1 ~ I < s and js is mapped into js — 1, from 



equation (3.12) the operator S always decreases the ranking of the terms on which it 
acts. 



This means that the residual term into the equation (3.11) only contains coeffi- 
cients in /3 (i.e /3) already determined when considering the corresponding equation 
at the higher orders. Hence we can solve the linear equation(3.11 ) w.r.t 



in terms of and in terms of already determined higher ranking coefficients. 

This concludes our proof. 
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4 Classification results 



We now proceed to the classification of involutive 1-forms 



, def 



/■'^ = /(«) + eAu^ + (fiiM^.^. + B2UI) + (CiM^^^^a. + C2U^u^^ + C^ul) 
+ {DiU4x + D2UxUxxx + Dsul^ + Diuluxx + D^ul) 

+ {Eiu^x + E2UxU4x + E^UxxUxxx + E4uluxxx + E^Uxul^ + EquIuxx + E-jul) 

+ ... 

(4.2) 



up to a certain order w.r.t. tlie Poisson bracket defined in (2.5), that is 



In virtue of the theorem (3.2) we can take in its normal form (4.1) without loss 
of generality. We have the following 

Theorem 4.1 If a{u) 7^ 0, then up to 0(e®), the normal forms of integrahle conser- 
vation laws 

Ut = dxUjt^^ 



and their commuting flows 



Ur = dxUjf^ 



are parametrized only by the functional parameter a{u) (here uj'^^J and u'j^^ are re 



spectively given by (4-1) and (4- -2)) 



Proof: The constraint that the Poisson bracket {u'^^t-^ , u'j-'^^ } vanishes up to the 
order O(e^) provides the following set of conditions on the deformation coefficients 

0(e°) — )■ no conditions 

0(e) — )■ Gives A in terms of a and /: 

A{u) = la{u)nu) (4.3) 
O(e^) —7- Gives Bi and B2 in terms of hi, f and a{u): 
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O(e^) — > Gives Ci, C2 and C3 in terms of the small letters and / 
C, = ia^a'r + -Ci/" + ^aV^'^ 

together with the constraint 

h = (aV2!)' (4.4) 

O(e^) Gives Di, D2, -D3, -D4 and 1)5 in terms of the small letters and / 

D^ = la'a'f^'^ + -a'a"r + —a^'^f^'^ + -a'{a'fr + -dj" 
8 6 120 2 ^ ^ 2 ■' 

-D2 = —a a f"- ' + -d2a f +-aaaf + dif + -a a f + J + a > + 
16 2 4 6 48 8 

Ds = ^a^a'a'T + ^a'a-f + ^«'«'7^^^ + ^«(«') V'" + ^a^(a')V^^^ + ^^1/'" + 

+ -a{afa"f"' + Aa^a'a"f^^ + dj^^^ + -a2(a') V^'^ + -a^a"'f^^ + 
10 8 3 



+ 

and the constraint 



^aiara"'r+la{a'ra"f^'\ 

ci = (aV3!)" (4.5) 
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O(e^) — )■ Gives Ei, E2, E3, E^^, E^, Eq and E^ in terms of the small letters and / 
(whose cumbersome expression is not presented here) and the constraints 

d, = (aV4!)'" 

24 3 ^ ^ ^ ^ 18 

■ 

We call integrable viscous conservations laws all conservation laws obtained from 
this deformation procedure extended at any order in the parameter e and such that 
a{u) 7^ 0. The perturbative calculations performed above support the formulation of 
the following 

Conjecture 4.2 The normal form of integrable viscous conservation laws associated 



with the 1-form (4-1) as well as their commuting flows associated with the 1-form {4-2) 
are uniquely determined at any order in e by the non-vanishing functional parameter 
a{u). 

It should also be noted that if a{u) is allowed to be vanishing the deformation proce- 
dure develops a branching. The branch a{u) = corresponds to the case of dispersive 
perturbations. In particular, Hamiltonian dispersive perturbations have been exten- 
sively studied in a number of papers (see e.g. [131 1121 [23]). It is conjectured that 
all deformations are uniquely specified by a number of functional parameters that 
in the bi-Hamiltonian setting are called central invariants [H]. The functional pa- 
rameter a{u) plays in the present context the same role as central invariants in the 
Hamiltonian setup. Hence, we refer to a{u) as viscous central invariant. 

An important well-known example of integrable viscous conservation law as been 
already mentioned above and it is given by the Burgers equation 

ut = dx {u^ + eux) = 2uux + eun,x, (4.6) 

Let us observe that Burgers' equation already appears in its normal form and it pos- 
sesses a constant viscous central invariant a{u) = 1. Conversely, one can immediately 
check that in the case of constant viscous central invariant the 1-form u'^t'^ truncates 
at the first order in e providing the Burgers equation. Higher flows of the Burgers 
hierarchy can be obtained as a deformation of higher flows of the Hopf hierarchy 



obtained choosing in (4.2) f{u) = m", n > 2. 



5 Linear viscous central invariants 

5.1 The viscous analogue of the Camassa-Holm equation 

In the present section we show that the deformation procedure discussed above pro- 
vides, in the case of the linear viscous central invariant a{u) = u, the following 
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non-evolutionary conservation law 

Ut — eUxt = dx {u^ — euUx) 



(5.1) 



Equation (5.1) first appeared as a scalar reduction of two-component Camassa-Holm 
type equations [311 [13 El] . 

This equation belongs to the class of viscous conservation laws under consideration 
and can be written in evolutionary form as a formal series in e. In fact, the application 
of the formal inverse operator 



(1 - ed^r' = 1 + ed, + e'dl + ... = J] 6^9 



(5.2) 



k=0 



to both sides of (5.1) gives 

ut = {I - edr,)~'^[2uuj: - e{uua,x + ul)] 



il-edx)-'d.. 



u 



+ (1 - tdx) - 



u 



dx 



+E^'3i(Y) 

fc=0 ^ ^ 



(5.3) 



We show that equation (5.1), or equivalently (5.3), can be brought to the normal 
form via a Miura transformation. This statement is made precise by the following 

Theorem 5.1 The equation (5.1) is reduced to its normal fo 



form 



via the Miura transformation 



(5.4) 
(5.5) 



where v = v{x,t) is a solution to the equation (5.1) or, equivalently, (5.3). 
Proof: Inverting the Miura transformation (5.5) we have 

V = {l- edx)'^u. 
Differentiating by t and using the equation (5.4) we get 



[I - edx) ^ut = (1 - edx) ^d^ j ^ e'^uu^k) ] = 9^(1 - ed^) ^uv 



dx{l - edx)~\v^ - evv^x) = dx ^'^l 

k=0 



3\1 



°° /2\1 r2 °° /2\ 

k=l ^ ^ J L A:=0 ^ ^ 



where we used the expansion formula (5.2). The theorem is proved. 
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Comparing with the results of the previous section we see that up to the fifth 
order in the deformation parameter the equation (5.4) coincides with the equation 



def 

2 



with a{u) = u. According with the Conjecture 4.2, the equation (5.4) is expected to 



be the only scalar integrable viscous conservation law, up to Miura transformations, 
possessing a linear viscous central invariant. We note that the positive hierarchy 
is obtained as the set of symmetries corresponding to the choice f{u) = u"^, n = 
0,1,2,3,... (including the original equation) while the negative hierarchy represents 
the set of symmetries corresponding to the choice f{u) = n = 1, 2, 3, .... In the 
next subsection we write explicitly the negative and positive hierarchies and provide 



an alternative proof of the integrability of equation (5.1) via linearization 



5.2 The negative hierarchy 

The analysis of deformations discussed in the Section |4] suggests that the normal form 
of the equation (5.1) is not truncated. Nevertheless, there exists a family of fiows in 
the hierarchy possessing a finite deformation, that is the negative hierarchy 

1' 



u 



1 



n 



1,2,3,.... 



The first equation of this hierarchy obtained for n 

ut_i = dx 



ed^- 



u 



u 



(5.6) 



(5.7) 



has been included by F. Calogero |1] within a list of nonlinear PDEs that are in- 
tegrable via a nonlinear change of variables. Subsequently, it was observed in [31] 



that the equation (5.7) can be mapped into Burgers' equation. Let us note that the 
hierarchy (5.6) can be written equivalently in terms of a recursion operator as 



where 



is the inverse of the operator 



dx{l 



ed^)-d,. 



-1 



u 



such that equation (5.4) can be written as 



Rux- 



{5.t 



(5.9) 



(5.10) 
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Operators (5.8) and (5.9) can be viewed viewed as a reduction of the recursion op- 
erators for, respectively, negative and positive hierarchies of the system studied in 
[in]. A direct proof that (5.8) and (5.9) are indeed recursion operators is based on 
the apphcation on a criterion proposed by P. Olver in [30] . 

We now prove the integrabihty of equation (5.1) in two steps. First, we show 
that the flows (5.6) can be mapped to the heat hierarchy by a sequence of reciprocal 
and nonlinear transformations and as a consequence commute. Secondly, we prove 
that the equation (5.4) transforms into the Klein-Gordon equation and consequently 
commutes with the flows (5.6). 



Lemma 5.2 The family of flows (5.6) is transformed into the Burgers hierarchy via 
the reciprocal transformation 



v{ip,t. 



J t—n) 



-d^x{(p,t_i, . . . ,t_ 



where 



u dx, 



(5.11) 



that is the function v is a simultaneous solution to the Burgers hierarchy 



n = 1,2,3, . 



Proof Introducing the potential if such that u = (fx, equations (5.6) after the 
integration w.r.t. x give 



edx- 



1 

For n = 1 the above equation takes the form 



(!)• 



1 , ^xx 
— +€ — . 

^x 



Introducing the reciprocal transformation of the form x = x{!f,ti,t2, 



Vx 



X 1/ 



Vx 



X 



X,, 



X// 



n = 1,2, 



(5.12) 

(5.13) 
one has 
(5.14) 



Using relations above into the equation (5.13) one arrives to 

Xf_-^ ( ~\~ €(9(^^ 

that is just the potential Burgers equation, i.e. v = —x^ satisfies the Burgers equation 

Vt^i = d^{v + ed^) V. 

Let us now proceed by induction assuming that the proposition is true for the 
n— th equation and consider the (n -|- 1)— th equation 

n+l 



1 



^x 



(1) 



1 

^x 



1 
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Using the relations (5.14) and observing that dx = x^^d^, we arrive to the equation 

xt.„_, = {x^ - ed^) xt_„. 

Since 



Xt. 



^ Xtp ~t~ adi^ 



we have 



xt_„_, = - {-x^ + ed, 



\n+l 



\p J X^p 



Finally, let us differentiate the expression above by if and substitute x^p = —v. The 
lemma is proved. ■ 



Theorem 5.3 The flows (5.6) commute 



Proof: In virtue of the Lemma (5.2) the set of flows (5.6) is transformed into 



Burgers' hierarchy via the reciprocal transformation (5.11). On the other hand it 
is known that the Cole-Hopf transformation v = ed^ log w brings Burgers' hierarchy 
into the heat hierarchy 

This result can be straightforwardly verified for n = 1 and then proved by induction 
for any n. Flows of the heat hierarchy clearly commute, i.e. dt^„dt_^w = dt_^dt_„w = 
rpj^^ theorem is proved. ■ 



The following theorem uncovers the relation existing between the equation (5.1) 



presented above and the integrable hierarchy (5.6). 



Theorem 5.4 The flow of conservation law (5.4) commutes with all flows of the 



hierarchy (5.6). 



Proof: Let us note that the equation (5.4) can be written in the more compact 



form 



vt = dx [v{l - edx) ^v] . 



(5.15) 



Introducing the variable (p defined by (fx = v, integration with respect to x provides 
us with the equation 



(Ptifx - (fl+e {ifit^xx - ^X^Xt) = 0, 



(5.16) 



where we used the fact that the integration constant can always be eliminated by 
a shift of the independent variable (p — )■ (f + f{t). The reciprocal transformation 
X = x{ip,t) brings the equation above into the form 



ex 



(ft '^ip 



X(nXf 1, 
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and the change of variable x = —e\ogw maps the equation above into the Klein- 
Gordon equation in light cone variables 



e w^t + w = 0. 



Note that both the hierarchy (5.6) and the equation (5.4) can be linearized via the 
same change of variables that brings them to the heat hierarchy and Klein-Gordon 
equation respectively. Since the Klein-Gordon equation is compatible with all mem- 
bers of the heat hierarchy, the theorem is proved. ■ 



5.3 The positive hierarchy 



Given the equation (5.4) in the form (5.10) it is natural to consider the family of flows 
of the form 

Ut„ = R^u,, n = l,2,... (5.17) 



where R is deflned in (5.9). We now prove that the family of flows (5.17) can be 
mapped into a family of linear commuting flows whose flrst member is the Klein- 
Gordon equation and for this reason we name it Klein-Gordon hierarchy. We have 
the following 



Theorem 5.5 The family of flows (5.17) is mapped to the Klein-Gordon hierarchy 
that is given by the following set of linear commuting flows 



n = l,2,... 



{5.U 



via the nonlinear change of variables 

x{ip,ti,t2, 

where 



-e log w 



X = (p 



and 



u dx. 



Moreover, the positive hierarchy (5.11) and the negative hierarchy (5.6) commute. 

Proof: Let us flrst note that the theorem holds true for n = 1 in virtue of the 
proof of the theorem (5.4). Let us now consider the equations (5.17) for n > 2 in the 
form 

Ut„ = Rut„-„ n = l,2, .... 

Introducing the function (f such that u = (px, following the same steps as in the n = 1 
case the above equation can be equivalently written as 
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Introducing the reciprocal transformation 



X = X{ip,ti,t2, ...) 

we arrive to the equation 

ex^t„ - x^xt„ + xt„_, = 0. (5.19) 
Setting X = — elogw one readily gets the following linear equations 

ew^t^ + wt„_^ =0, = 2, 3, . . . 



that can be equivalently written in the evolutionary non-local form (5.18). It is 
immediate to check the commutativity condition for the flows of both the positive 
and the negative hierarchy, i.e. dt„dt,^w = dt^dt„w for any n,m & Z. The theorem 
in proved. ■ 
Summarizing, the normal form of equations associated with a linear viscous central 
invariant can be mapped into the Burgers equation by means of a reciprocal transfor- 
mation and then to the heat equation via a Cole-Hopf transformation. Interestingly, 
the same sequence of transformations brings on one hand the negative hierarchy into 
the positive Burgers hierarchy and then to the heat hierarchy. On the other hand, 
the positive hierarchy is mapped into the negative Burgers hierarchy whose poten- 



tial form is given by the recursive nonlocal equations (5.19) that is linearised into 
the Klein-Gordon hierarchy (see also Figure 1)). Hence, the Klein-Gordon hierarchy 
can be interpreted as the negative heat hierarchy. Indeed, introducing the recursion 
operator 

7^ = ed^ 

the Burgers' hierarchy is written as 

= n^'w^ n = 0, 1, 2, . . . 
while the Klein-Gordon hierarchy is 

w. = n-'^w n = 1, 2, . . . . 



6 Local well-posedness of the Cauchy problem 
6.1 Kato's method 

The present Section is devoted to the study of the local (in time) well-posedness of 



the periodic Cauchy problem for (5.1) by using Kato's method [18] that is here briefly 
outlined. 
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Miura transformation 


= dx (Xlfc ^^v,U(^k)) ^ Normal form 


Negative hierarchy 
i 


Positive hierarchy 
i 


i 


i 


Heat hierarchy 


Klein-Gordon hierarchy 



Figure 1: Positive and negative integrable hierarchies associated with viscous con- 
servation laws possessing a linear viscous central invariant. 



Let us consider the Cauchy problem for the quasi-linear equation of evolution: 

+ A{v)v = f{v), t>0 t;(O)=0. (6.1) 

Let X and Y be reflexive Banach spaces with Y continuously and densely embedded 
in X, let 5* : y — X be an isomorphism and suppose that the norm of Y is chosen 
in such a way that 5* is an isometry. Assume now that the terms appearing in (6.1) 
satisfy the following conditions: 

1. A is an operator-valued function defined on Y, that is for each y & Y, A{y) : 
D{A) (Z X ^ X is a linear operator on X (in the interesting cases A{y) is 
an unbounded operator). Assume that A{y) is quasi-m-accretive, uniformly for 
y & Y with \\y\\Y < M. In other words, for every constant M > 0, there is a 
real number u such that for every y ^ Y, with ||y||y < M, —A{y) generates a 
C°-semigroup {e~^^^'^^t>o with 

||e-*^(^)|U(x) <e-*, t>0, 

where || ■ \\c{x) is the operator norm. 

2. For each y &Y, A{y) is a bounded linear operator from F to X and moreover 

\\{A{y) - A{z))w\\x < fJ'AWy - z\\x\\w\\y, y,z,weY, (6.2) 
for some constant fiA depending only on max{||y||y, ||-2||y}. 

3. for any M > 0, the inequality 

\\{SA{y) - A{y)S)S-'w\\x < fiiiM)\\w\\x, yeY, \\y\\y < M (6.3) 
holds for all w eY, where fii{M) > is a constant. 
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4. For each M > 0, / is a bounded function from {y : \\y\\Y < M} to Y. Also 
we have 

\\f{y)-f{z)\\x<fi2\\y-z\\x, y,zeY, (6.4) 

and 

\\f{y)-f{z)\\Y<f^3\\y-z\\Y, y,zeY, (6.5) 

for some constants fi2 and /is, where /i2 depends only on max{||y||2, ||-2||z} and 
II3 depends only on max{||?/||Y, where Z is a space such that Y C Z C X , 

all with continuous inclusions. (In our application it will be X := L^,Y : = 
H^, Z := L°° all for functions defined on the circle.)]^ 

Then one has the following Theorem (see (18]): 

Theorem 6.1 Assume conditions (1), (2), (3) and (4) above hold. Then for any 



(f) E Y , there is aT > 0, depending only on ||0||y and a unique solution v to (6.1) 
such that V G C°([0,T],F) fl C^{[0,T], X) . Moreover, v{t) depends continuously on 
the initial data = t'(O) in the Y-norm. 



6.2 The periodic Cauchy problem for the equation (5.1) 



Let us consider the equation ( |5.1[ ) with fixed e and, for convenience, we set e = 1 and 



replace t — ?■ —t. With these notations the initial value problem for the equation (5.1) 
will read as follows 

rrit = —vrrix — v^m v{x, 0) = w{x) (6.6) 
where m := v — Vx and w : §^ — )■ M. The above time inversion will result into proving 



the local well-posedness of the Cauchy problem for equation (5.1 ) backward in time. 



We will finally observe that the result holds also true for the forward evolution. 



Using the fact that v = S where S = 1 — dx the problem (6.6) can be 
equivalently formulated as follows 

rrit = —{S^'^m)mx — m{S~^m)x, m(x, O) = 0(x) (6.T) 

where : §^ — )■ M. Verifying the conditions for the validity of Kato's theorem will 
guarantee well-posedness for G H^(S^). 

From now on all the functional spaces introduced are referred to functions on the 
circle and we for sake of simplicity will be denoted just as L^, H^. 

We have the following 



^Note that these assumptions are exactly hke those stated in [7] , except for the Lipschitz inequahty 



(6.4). Indeed in it is required that ^2 depends only on niax{||y||x, H^ljx}, while here we require 



dependence on a bigger quantity. As it can be seen in [18j page 40, second inequality from above, 



this is not affecting the statement of Theorem 6.1 in our case. 
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Theorem 6.2 Let G H^, then there is aT > 0, depending only on \\4>\\h'^ such that 
(6.7) has a unique solution m G C°([0, T], /f-*^) fl C"'^([0, T], L^) and moreover m as a 
function oft depends continuously on cf) in the H^-norm. 



Proof: We have just to check assumptions (1) to (4) of Theorem 6.1, We choose 
as X = and as y = H^. On we use the usual L^-norm. Recall that is 
continuously and densely embedded in L^, since for any y G WvWl'^ ^ WvWh^ (by 
definition \\y\Wi ■= \\y\\'^2 + WyxWi^)- Also we choose as an isometric isomorphism 
S = 1 — dx, so that 

\\Sy\\x = \\y-yx\\h =<y -yx,y -Vx >l^= 

= \\y\\h + Ib:r||i2- < yx,y >L2 - <y,yx >l2, 

where the terms < y^, y >l2 + < y,yx >l'^ give a zero contribution due to integration 
by parts which holds for functions in (see for instance Corollary 8.1 [3]). Therefore, 
with the choice of S as an isometric isomorphism, the norm on is equal to the 
standard norm on H^. Comparing (6.7) with ( 6.1[ ), we choose 

A{y) = {S-^y)dx, fiy) = -yiS-'y)x, yeY = H\ 

We start verifying assumption (1), fixing M > and y eY with ||?/||y < M . First 
notice that {S~^y) G Y but since 

iS-'y)x = S-'y - y (6.8) 

it turns out that S^^y has actually a representative (this is because is em- 
bedded in C°). In order to check condition (1) then, we need just to show that 
suPa;e[o,i] \ i'^~^y)x\ is bounded by a constant depending only on M (see [18J, p. 38). 
On the other hand, by (6.8) 

sup \{S~'^y)x\ < \\S~^y\\Lo^ + Wyh-- < c\\S~'^y\\Hi +c||?/||^^i, 
xe[o,i] 

where c is the constant entering in the Sobolev inequality corresponding to the inclu- 
sion )■ L°°. Since S is an isometry we have ||S'~"^?/||iifi = ||?/||l2 and ||?/||l2 < WvWm 
from which we deduce 

sup \iS-'^y)x\ < 2cM. (6.9) 
ie[o,i] 

Therefore, A is quasi-m-accretive on with u bounded by cM. 

To check assumption (2) we first show that A{y) is a bounded operator from 
to for every y G H^. Indeed we have, for w G H^: 

\\A{y)w\\L2 = \\{S''^y)wx\\L2 < Ik^^lUa < c\\S'^y\\Hi\\w\\m = 
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To conclude the verification of (2), namely showing that (6.2) is fulfilled, we just 
observe that 

\\{A{y) - A{z))w\\l^ = \\S-\y-z)w4L^ < \\S-\y - z)\\l^\\w^Al^ < 



< c\\S {y - 2;)||/^i||u;||j/i = c\\y - z\\l2\\w\ 



Let us now verify the estimate (6.3). We notice that 

SA{y)v - A{y)Sv = -{S'^y)^v^, 
then, for all y,w & H^, with < M we have 

WiSAiy) - Aiy)S)S-'w\\L^ = \\ - iS-'y).,iS-'wyAL^, 



which in virtue of (6.8) becomes 

\\{S-'y-y){S-'w-w)h2. 

Therefore, 

WiSAiy) - Aiy)S)S-'w\\L^ = \\iS-'y)iS-'w) + yw ~ iS-'w)y - iS-'y)wh^ < 

< \\iS-^w)iS-'y-y)\\L2+\\wiS-'y~y)\\L^. < {WS-'wh^ + Mr^) {\\S-'y\\L^ + Ml^) . 
Since HS^^ylU- < c\\S-^y\\m = c||?/||l2 < c\\y\\m, llvh^ < c\\y\\m and ||S'"^w||l2 < 



ll-S" = IblUa, we get 

WiSAiy) -Aiy)S)S-'w\\L^<Ac\\y\\Hi\\w\\L2<AcM\\w\\L^. 

In order to check the assumption (4), we first show fiy) := —yiS^^y)^ is bounded 
from L:= {y eH^ : WvWh^ < M} to H\ We have 



US-'y)xWm = \\yiS-'y)xWl^ + {yiS-'y)^) || 



2 

L2) 



and 



2 11^4 



\\yiS-'y)Al^<\\iS-'y)A\U\y\\l^<'^c^M' 



in virtue of (6.9) and ||?/||l2 < ||?/||_f/i < M. To control the other term, we observe 
that 

Wd. (yiS-'y),) IU2 < Wy.iS-'y).WL^ + WyiS-'y).M^ < 
< WiS-'y)M-WyM'^HyiS-'y-y-y,)WL^<2cM^+WyWL^iWS-'yWL^H^^^ < 

< 2cM^ + cM(3M) = ScM^, 
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where we have used once again (6.9) and ||?/a;||L2 < WvWh'^ < Combining these 
estimates we obtain 



Let us prove the first Lipschitz inequahty (6.4) is fulfilled. We have 



\\f{y) - f{z)\\L^ < \\yS-'y - zS~'z\\l2 + \\y' - 



L2. 



For the first term \\yS ^y — zS ^z\\l2, adding and subtracting zS ^y we obtain 

\\yS~^y - zS~'^z\\l2 < WS^'^yWio^Wy - z\\l2 + \\S'^{y - z)\\l^\\z\\l2 < (6.10) 

< ||S'"^y||//i||?/ - z\\l2 + \\S'^{y - 2;) II //I II 2; 11^2 = (||y||L2 + ||2;||i2)||y - 2||j:2 (6.11) 

< fxWy - z\\l2, (6.12) 

where /i = max{||?/|| ^^2, || 2; 11^,2}. Notice now that we can not bound the term Hy^—z^ 11^,2 
with a Lipschitz constant depending only on the norms \\y\\L^ and ||;2||l2, but we will 
have to use the L°°-norms that are bigger. This is still sufficient, as one can see in 
the second inequality at page 40 of [18]. Then it is immediate to get 



z IU2 < (||y|U°° + lkllL°°)||y - 2;|U2- 



(6.13) 



Combining this with the inequality controlling \\yS ^y — zS ^z\\l2, (6.4) is checked. 
Finally we prove the last Lipschitz inequality (6.5). We have for every y,z & 

\\y{S-'y), - z{S-'zUhi < MS-'y) - ziS-'z)\\H^ + \\y' - z'\\h^, (6.14) 



using (6.8). We first control the first term on the right: 

MS-'y) - ^(5-^^)11^^! < US-'y) - ziS-'z)\\L2 + \\d, {y{S-'y) - ziS^'z)) 1^2 < 

< MS-'y) - z{S-'z)\\l2 + \\y.iS-'y) - z,iS-h)\\L2 + MS-'y), - ziS-'zU\L2 < 
< + |k||L2)||?/ - z\\l2 + \\y^{S'^y) - z^{S''^z)\\l2 + 

+i\\y\\L2 + \\z\\L2)\\y - z\\l2 + {\\y\\L-o + ||2;|U->)||y - 2;||l2, 



using the inequalities (6.10), (6.11 ) and (6.13). Therefore we obtain for the first term 
on the right hand side of (6.14) 



US-'y) - z{S-'z)\\Hr < 3i\\y\\L^ + \\z\\Lo.)\\y - z^r + \\y.{S~'y) - z,{S-'z) 



|L2. 



On the other hand, adding and subtracting yx{S ^z) inside ||2/x(5' ^y) — z^^S ""^ 2;) 11^2 
we obtain 



\y,{S-^y)-z,{S-^z)\\L2 < \\y,{S-\y - zmL2 + \\{S-'z){y, - z,) 



< 
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< IIS ^{y - z)\\L^\\yx\\L2 + \\S ^z\\Loo\\y^ - Zr,\\L2 < 
< c\\S''^{y - 2;)||j^i||?/||//i + c||S'"^2||/^i||y - z\\hi < c(||?/||j/i + WzWh^) \\y - z\\hi. 



This proves that the first term on the right hand side of (6.14) is bounded by n\\y 



zWh^ where /i depends only on max{||?/||/^'i, ||2;||//i}. Finally to control the second 



term on the right hand side of (6.14) observe that 



lb' - z^Wm <\\iy + z){y - z)\\l2 + \\d, {{y - z){y + z)) < 
<\\{y + z){y - z)\\l2 + \\{y- z),{y + z)\\l2 + \\{y- z){y + z),\\l2 < 

< \\y + ^llL°°l|y - ^IIl2 + ll(y - 2:)x||L2||y + ^IU°° + \\y - ^IU°°ll(y + < 

< c(||y||/fi + |k||//i)||y-2;||Hi+c||y-2;||j/i(||y||/fi + ||2;||//i)+c||y-2;||//i(||y||^^i + |^ < 

< 3c(||y||^^i + ||5;||^i)||?/ - 2;||^i. 
Therefore we can conclude that 

\\y{S'^y)^ - z{S''^z)^\\hi < yUsllz/ - ^ll^i, 

where /is depends only on max{||?/||j:/i, ll-^H/fi} so assumption (4) is satisfied and the 
Theorem is proved. ■ 



Remark 6.3 In general the time T in the Theorem 6.2 depends on the norm \\(j)\\H\, 
hence for each choice of M > such that \\4>\\h^ ^ M we will have a corresponding T 
depending on M (see fTE^)- Moreover, since the operator A{y) used in the above proof 
generates a C'^ -semigroup, and not just —A{y) (see for instance f3^ . Theorem 12.26) 



we have local well-posedness also forward in time for the original equation (5.1) 



7 Numerical solutions 

The present section is devoted to the numerical study of the initial value problem for 



the equation (5.1) in the case of periodic boundary conditions. Let us observe that 



introducing the auxiliary variable P{x,t) such that 

(1 - a.)p = lv\ 



equation (5.1) can be equivalently written as follows 

(1 - dx) {vt - vvx - Px) = 



(7.1) 



(7.2) 
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Since the only smooth function / fulfilhng periodic boundary conditions such that 
(1 — dx)f = is f{x, t) = 0, the equation (7.2) is equivalent (on smooth solutions) to 
the system 



R =P - ^v^ 



Vt =VVx + Px, 

1 . (7.3) 
2 

Numerical solutions are obtained by using a pseudo spectral method implemented in 
the NDSolve Package of Mathematica 9 [34J. Although local existence of solutions 
is guaranteed by the Theorem 6.2\ we observe that positive definite initial data seems 
to support global existence (in time) of solutions. Contrarily, finite-time break up 
occurs in other cases. In Figure [2] we compare the evolution at different time steps of 
the positive definite sinusoidal initial data 

vi{x) ■.=v{x,0) = sm[ — ) +2 (7.4) 



12/ 

/ Tlx \ 

V2{x) ■.=v{x, 0) = sin i^—J + 2. (7.5) 

From the definition of P{x, t) we have respectively 

, ^, , 9 9cos(|7rx) 37rsin(7ra;) 247r cos (p^vrx) 288 sin (^p^vrx) 

P, r) ■=P r 0) = — -\ — - -\ — -\ — 

r^yx). r^yx,u) ^ ^ ^ 2{36 + tt^) 144 + tt^ ^ 144 + vr^ 

, , 127rcos(^) 9cos(f^) 72 sin (^) 37rsin(^) 9 
r^2[-L) . ryx,yi) + gg 4(7^2 + 9)+ 7^2 + 35 +9(^^2 + 9) +4- 

In both cases, after the initial steepening both solutions tend to vanish in time. It 
should be noted that the steepening is more pronounced for the datum of larger 
period (7.4) as expected, although the effect of dissipation prevents it to reach the 
gradient catastrophe. In the case of the Camassa-Holm equation it was proved that 
if the quantity u{x, 0) — u^^(x, 0) does not change sign, the periodic Cauchy problem 
is well-posed in |[7J. It would be interesting to prove a similar results for the 
equation (5.1). 

A completely different scenario takes place when considering initial data that 
change sign, as for instance an odd initial datum as 

v^{x) := v{x, 0) = sin j . 

The corresponding initial function P is 

^ ^ 37rsin(^) 9cos(^) 1 

As Figure [3] shows no damping occurs in this case and the solution seems to develop 
a gradient catastrophe in finite time at the inflection point a; = 0. A similar phe- 
nomenon has been observed in the case of Camassa-Holm equation for which there 
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X 



X 



t = 5 



Figure 2: Time evolution for tlie initial datum v{x, 0) = sin (ff) + 2 (solid line) and 
0) = sin + 2 (dashed line) with periodic boundary conditions. 



v[x 



exists a class of initial conditions (m(x, 0) G H'^ with u{x, 0) odd and u{x, 0) ^ 0) 
such that the solution blows up in finite time[6l [7]. 

Remark 7.1 Let us observe that ifG is the Green function for the operator (l—dx) on 



the circle, then from the definition of the funtion P (7.1 ) we have P = G-k ^v^, where 



•k denotes the convolution product. In full analogy with the C amass a- Holm equation, 



using this representation for P the equation (5.1) can be written as a conservation 
law with nonlocal flux 



2 ' V 
2 



8 Quasitriviality, transport equations and deformed 
hodograph formula 

A Miura transformation is defined as a change of variables of the form 

u{v,v^,...) = ^e''Fk{v,v^,v^^,...), degFk = k (8.1) 

k 

where coefficients Fk are differential polynomials in the derivatives Vx, v^x, ■■■ A natu- 



ral generalisation of (8.1), referred to as quasi-Miura transformation, is obtained just 



relaxing the assumption on the polynomial form of F^. 
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Figure 3: Steepening of tlie solution for tlie initial datum v{x, 0) = sin (y|) . 



According to the results of [23j, any evolutionary PDEs of the form 
vt = vvx + ^ t^Pk{v, fx, fxx, ■■■), degPfc = k 

k 



can be reduced by a quasi-Miura transformation to the Hopf equation (3.2 ). Moreover 
any evolutionary symmetry of the equation 



k 

is reduced by the same transformation to its dispersionless limit 

Ur = c{u)Ux- 

In particular, given an integrable hierarchy of evolutionary PDEs, there exists an 
invertible quasi-Miura transformation that brings any equation of the hierarchy to 
the corresponding equation of the Hopf hierarchy. 

The method for constructing recursively the terms of a quasi-Miura transformation 
at any order in e has been illustrated by S. Liu and Y. Zhang in [23] and it is based on 
the construction of infinitesimal generators Xi, X2, ... for the reducing transformation 

v{u) = expX{u) = u + X{u) + ^X (^X(u)) + (8.2) 
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where 



and 



OU OUr 



X = eXi + 6^X2 + .. 



In the following, we propose an alternative but equivalent version of the Liu-Zhang 
approach that is based on the direct solution of transport equations for the class of 
equations of the form 



Looking for asymptotic solutions in power series of e of the form 



V = u + ev'^ + e^v^ 



1.3) 



where u = the equation (8.3) splits into a quasilinear equation for uq plus a set of 
transport equations 



Lu = 



1.5) 



where 



L = dt — 2udx L* = dt — d^lu 



are, respectively, the first order differential operator of the Hopf flow and its adjoint. 
Hence, the problem of constructing the reducing quasi-Miura transformation is re- 



duced to the classical problem of constructing the expansion of the form (8.4) via the 



transport equations (8.5). Such solutions are referred to as formal solutions to the 



equation (8.3). Let us observe that once a solution u to the first equation in (8.5) 
is known, higher order corrections f ^, f^, ... are obtained by solving a sequence of 
linear PDEs with coefficients depending on u. Hence, the asymptotic formal solution 



of the form (8.4) can be interpreted as a transformation that brings any solution to 
the Hopf equation to a solution of the deformed equation (8.3) [13] . It turns out that 



for the class of equations under consideration such a transformation is of quasi-Miura 
type. Given a solution to the Hopf equation Lu = that is implicitly given in terms 
of hodograph equation: 

x + 2ut + f{u) = 0, (8.6) 
the n-th transport equation can be written in the form 
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where we have observed that using the differential consequences of the equation (8.6) 
u.. = f"{u)ul u^,^ = f"\u)ul + ?>{f"ful etc. (8.8) 
the r.h.s can be written as a suitable function of u and Ux only. The general integral 



of the equation (8.7) is readly obtained as a function of the variables u and of the 
form 

V'^ = Pniu,Ux)Ux + hn{u,Ux) (8.9) 

where 

P„(m, Ux) = J —Fn{u, (p) dip 

and hn{u,Ux) is the general solution to the homogenous linear equation L*hn = 0. 
Using the method of characteristics one can show that 

hniu.Ux) = gn{u)Ux 

where gn{u) is an arbitrary function of its argument. 

For the Burgers equation, that is obtained from (8.3) by choosing a constant 
central invariant a{u) = 1 the set of transport equation takes the a simple ricursive 
form 

L*v'^ = dx{ux) 

= dx v'v"^'' + j , iV = 2, 3, . . . 

with the notation = u. 

Using the argument above one can prove that the solution to the ri— th transport 
equation for the Burgers equation is written as a polynomial in Ux of the form 

2n-l 

t;" = ^ an,n+j + gn{u)ux (8.10) 
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where coefficients a are determined by recursion as follows 

1 



_^ ["n-l,n+l + + ^n-l,l] 

+ (n + J - 3)(n + J - 1)(/"(m))2«„_i,„+,_3] , j = 3, . . . , 2n - 3 

an,3n-2 =^(^^"7^ [An-l,2n-2 + (3^ " 3)/"(M)A„_i,2n-3 + f^n-l,2n-3 

+ (3r2 - 5)(3n - 3)(r(7i))2a„_i,3,_5] 

lyin — 2) ^ ' -' 



where 



«i,2 = ^r(«) 



n-l 2i-l 



i=l j=l 

= (2n + 2j - 1)/" + (n + j - l)r'(^i)a„-i,„+,-i. 

Let us note that coefficients a's depend only on the variable u through the function 
/"(m) and its higher derivatives. The quasi-Miura transformation for Burgers' equa- 



tion can be recovered from the formula (8.10) by choosing gn{u) = and eliminating 



the dependence on function f"{u) and its higher derivatives via the triangular system 



of the form (8.8). It turns out that terms of quasi-Miura transformation at any order 
are rational functions of Ux, u^x etc. with no explicit dependence on u. For instance, 
up to O(e^), the quasi-Miura transformation for Burgers' equation takes the form 

I "^xx 2 f '^xxx '^xx\ I 3 f '^xx'^Ax "^xxx , 3'^xx'^xxx\ _i_ 4\ 

"""^ ^xJx '^^^Jx 

(8.11) 

A direct comparison with the Liu- Zhang approach to the solution of transport equa- 
tions [23 is readily made by observing that the functions f"{u), f"'{u), . . . play the 



role of parameters in the quadrature formula (8.9) as well as the functions Xuu = 
= f"'{u), ... do in Liu-Zhang's framework. 

Theorem 8.1 Given any solution to the Hopf equation 

Ut = 2uux, 
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via the hodograph formula (8.6), the formal solution of the Burgers equation of the 
form (8.11) satisfies the deformed hodograph equation 

x + 2ut + ujf = {) (8.12) 

where Uf is the deformed 1-form corresponding to a{u) = 1. 

Proof Let us consider the vector field 

Paf = d^{x + 2ut + f) = l + 2u^ t + dj\ (8.13) 

It can be viewed as a linear t-dependent combination of the vector fields 

d d d 

OU OU OU 



Applying to such vector fields the inverse of the quasi-Miura transformation (8.11) 
and taking into account the transformation law for vector fields 

x(u) ^ x(v) = + —a, + —ai + ...j A-(.),_,„,„.,.., 

we can immediately prove that 

• — )■ ^ as a consequence of the fact that |^ = 1. 



OU av 



• (<9x/) — (dxCOf)^- In the general case we only know that such a deformation 
exists up to the order e^, but in the case of Burgers it is defined for an arbitrary 
(analytic) function/. 

Combining the above results we obtain 

{l + 2uxt + dj)— ^ {l + 2vxt + dxUf) — 



This means that, given any hodograph solution u{x,t) the series (8.14) satisfies the 
equation 

l + 2vxt + dxUJf = 0. (8.15) 
Integrating with respect to x we obtain 

x + 2vt + Uf = c. (8.16) 

where c is a constant. Taking the limit e — it is immediate to check that the 
constant c must vanish. 
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Let us now consider the general case where the function a is not constant and 
suppose, as we have conjectured, that there exists an integrable hierarchy for any 
choice of the central invariant a{u). As a consequence of the results of S. Liu and 
Y. Zhang there should exist also in this case a reducing quasi-Miura transformation. 
One can easily check that such a transformation depends in general on v at any order 
in the deformation parameter. This immediately implies that the vector field is 
no longer invariant and an additional contribution must be taken into account. If 
this additional term is a total x-derivative, one obtains a correction to the deformed 



hodograph formula (8.12). Indeed if 

dv 
du 



l + d,{F{u)) 



then — > [1 + (F(-u(f )))] and therefore 

(l + 2M,t + 9,./)— ^ {l + 2v^t + d^{ujf + F{u{v)))) — . 

Proceeding exactly as in the above theorem one obtains the deformed hodograph 
formula 

X + 2vt + ujf + F = {]. (8.17) 
For example, in the case a{u) = u the reducing transformation is 



1 

u+-e 



1 u^u 



ui 



(4) ^ ^UUxxx 



M.7 



3 



7 uuz 



3 



+ 



uu 



2 



xxxIuMx MM^^ln-Ua- \ _|_ ^^^3^ 



whose inverse is 



u 



+ Vx In Vx 



1 v^v 



(4) 



1 



+ 



5 v'^v^ 



12 



T Vxx\y^Vx) Wa;a;mWa; 7 h" 7, 



1 v<x 
12 vl 



After some computations one gets 
'1 



dv 
du 

l + dx 



l + dx 
1 



-ehiUx + e' 



4 ul 4 Ur 2 u^. 



1 uujx 
3 ul 



-eln — 



12 t;3 



and therefore the correction to the deformed hodograph formula up to the secon order 
is given by 



F{v) 



1 



-e Inf^. 



1 V<a 

12 vl 
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The asymptotic approach based of the construction of formal solutions via quasi- 
Miura transformations or equivalently via transport equations is generally expected 
to provide an accurate local asymptotic description of solutions to equations of the 



form (8.3) for sufficiently regular initial data and sufficiently small times, i.e. before 
the generic solution to the Hopf equation Lu = develops a gradient catastrophe. The 
rigorous justification of this method requires a dedicated analysis that goes beyond 



the scope of the present work. However, for some equations of the form (8.3) rigorous 



results are already available (see e.g. As it was pointed out in [26J the fact 



that the quasi-Miura transformation (8.1) does not preserve the initial datum makes 
difficult to perform a direct comparison between the perturbed and the unpertuberd 
solution to a specific initial value problem. This problem can be readily fixed using 
the fact that terms at any order in e of a quasi-Miura transformation are particular 



solutions to transport equations (8.5 ) whose general integral is defined up to the kernel 
of the adjoint operator L*. Hence, given a solution to the Hopf equation Lu = with 
initial datum u{x, 0) = uq{x), the required perturbed solution v = u + ev^ + e^v^ + . . . 
that satisfies the same initial datum has to be such that v^{x, 0) = 0, f ^(x, 0) = etc. 
Such solutions are obtained by choosing the function hn{u,Ux) = gn{u)Ux in (8.9) in 
such a way that 



gniuo) = -pniuo,Uox) = -Pn Uq, -—, r 

/ (^^o) 



where we have used the formula 

1 



2t + f'iu) 



t=Q 



In the KdV case the existence of a transformation reducing KdV equation to Hopf 
equation and preserving the initial data was proved in [25]. The above arguments 
shows that, in general, the existence of this transformation relies on the freedom in 



the choice of solutions of transport equations (8.7). 



9 Dubrovin's Universality 

Let us consider the Hopf equation, that is the conservation law associated with the 
undeformed 1-form = 

ut = 2uux, (9.1) 
and its solution given in terms of the formula 

x + 2ut- f{u) = (9.2) 
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where the function f{u) is an arbitrary function that parametrizes the family of 
commuting flows associated with the unperturbed 1-form uj = f{u) 

Ut = f\u)ux- 

The deformation procedure proposed leads to the pair of involutive 1-forms 



def 



+ eaUx + e^aa'uxx + Oie 



UJ 



= /(«) + ^af'ux + e 



1 1 



/ £'// I ^ 'Z rill \ I ( / till I 2 ^(4) 

aaj + -a f ] u^x + { 7;(^a f + T^a ' 



6 



where a = a{u). Let us now proceed with the study of the critical behaviour of 
Burgers equation in full analogy with the case of dispersive Hamiltonian equations 
considered by Dubrovin in [12] . Introducing the deformed hodograph equation 



x + 2ut- cof^ = 



(9.3) 



where co'j:^^ has been specified for a constant central invariant a(u) = ao, we shall per- 
form a multiscale analysis about the generic point of gradient catastrophe {xo,tQ,UQ) 
such that 

xo + 2uoto-f{uo) = 2to-/'M=0 /"M = /'"M > 0. (9.4) 

Introducing displacement variables (x, t, u) as follows 

X = Xq + \°'x t = to + \H u = Uq + Xu (9.5) 

where A = e'^ is a small parameter, let us now expand in Taylor series the l.h.s of 
equation (9.3). Using the conditions (9.4) one obtains 



'"tttS 



6 



A+g(2-a) 



£111 



^2+g(l-2a) 



2 - 6 

= O (e4<?) + o (ei+«(2-.)^) ^ Q (^2+5(2-2.) ^) 

where the variable x is defined as 

X'x = + 2uoX^i = x-xo + 2uo{t - to), 



(9.6) 



with the notation /o = f{uo) etc. The request that all terms into the l.h.s of (9.6) 
contribute to the same order in e gives 



a 



Q 



Hence, we we end up with the expression 

~ I o — T ^ rill — 3 rill 

x + 2ut- -/o u - -ao/o u Ui 
6 / 



(9.7) 
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Then, the solution to the deformed hodograph equation in the vicinity of the critical 
point satisfies the second order ODE obtained by taking the limit e —t- in (9.7). The 
above equation can conveniently be written in non-dimensional form via the following 
rescaling of the dependent and independent variables 



X = siX t = S2T u = S3U 



where 



Sl 



a^f" 



1/4 



S2 



IsJ 
2 do 



1 

24 



1/2 



S3 



6ao 

f/r/ 

Jo 



1/4 



and U = U (X, T) satisfies the equation 

Uxx + SUUx + -UT = X. 



In terms of the original variables, the critical behaviour near the point of gradient 
catastrophe given a solution to the equation (9.8) is 

6 



U = Uq + 



6ao 

//// 
/o 



1/4 



3 fin 




"oi( 



X - xo + 2uo{t - to) 



p3/4 



24 



1/2 



(9.9) 

We also observe that the equation (9.8) can be transformed via a Cole-Hopf transfor- 

(9.10) 



mation of the form U = dx log w into the linear ODE 

wxxx - Twx = Xw. 



Equation (9.10) has the following general solution 

w{X,T) = ci oF2([l/2,3/4], ^ {T + xf)+C2 oF2([3/4, 5/4], ^ {T + Xf) (T + X) 

o4 o4 

+C3 oF2([5/4, 3/2], ^ (T + Xf) (T + Xf , 

where ci, C2, C3 are arbitrary constant of integration and 0F2 is a generalized hyperge- 
ometric function defined by 



+00 



oF2([a,/3], 2;) ■=Y, 



n=0 



'a)n{P)n n\ 



(9.11) 



where 



{a)n ■= (a)(a + l)(a + 2) . . . (a + n - 1), (a)o := 1 

and analogously for (/3)„. The function oF2([a, /3], z) defines an entire function when- 
ever a > and /3 > (it is clear from the expansion (9.11) that is dominated by the 
series of the exponential function up to a constant). 
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According with the conjecture formulated in [9j that generahses a result by A. 
Il'in [21], the critical behaviour (9.9) is provided by the particular solution to the 

Px{X,T) 



equation (9.8) 



U = dx\ogP{X,T) 



P{X,T) 



where P{X,T) is the Pearcey integral 



P{X,T) 



By a direct calculation using the identity 

/CO 
(-16^3 + ATz - 2X) e-(4^'-2T.2+2X.) ^ q 
-oo 



(9.12) 



one can directly verify that the function P(X,T) satisfies the linear ODE (9.10). 
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